Critical behavior of magnetic systems with extended impurities in general dimensions 
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We investigate the critical properties of d-dimensional magnetic systems with quenched extended 
defects, correlated in Sd dimensions (which can be considered as the dimensionality of the defects) 
and randomly distributed in the remaining d — Ed dimensions; both in the case of fixed dimension 
d — 3 and when the space dimension continuously changes from the lower critical dimension to 
the upper one. The renormalization group calculations are performed in the minimal subtraction 
scheme. We analyze the two-loop renormalization group functions for different fixed values of the 
parameters d,Ed- To this end, we apply the Chisholm-Borel resummation technique and report the 
numerical values of the critical exponents for the universality class of this system. 
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I. INTRODUCTION 

The study of critical properties of the various kinds 
of disordered systems remains one of the central prob- 
lems in condensed matter physics. We focus our atten- 
tion here to the case of e?-dimensional magnetic systems 
with quenched nonmagnetic impurities. The first ques- 
tion that arises here is: how does quenched disorder influ- 
ence criticality, i.e., will the universal scaling laws, gov- 
erning the critical behavior of the "pure" system be al- 
tered by quenched disorder introduced into the system in 
a quenched manner? 

The effect of weak quenched uncorrelated point-like 
disorder on the critical behavior of magnetic systems is 
predicted by the Harris criterion^: disorder changes the 
critical exponents only if the critical exponent a p of the 
pure (undiluted) system is positive: 



= 2 



dv p > 0, 



(1) 



v p being the correlation length critical exponent of the 
pure system. Of the d-dimensional spin systems that are 
described by the m-vector model, only the pure Ising 
model is characterized by a value of a p > and thus is 
affected by point-like weak disorder at criticality. 

Real magnetic crystals often contain defects in 
the form of linear dislocations, planar grain bound- 
aries, 3-dimensional cavities or regions of different 
phases, embedded in the matrix of the original crys- 
tal, as well as various complexes (clusters) of point- 
like non-magnetic impurities^. Systems with such 
"extended" (macroscopic) defects have attracted much 
interest 3 . 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 . 11 ! 12 ! 13 ! 14 . 

DorogovtseA* 3 - proposed the model of a d-dimensional 
m-component spin system with quenched random non- 
magnetic impurities, that are strongly correlated in e e i 



dimensions and randomly distributed over the remaining 
d — Sd dimensions. Such a system is no longer isotropic; 
the idea of two different correlation lengths naturally 
arises since the system is expected to behave differently 
along the directions "parallel" to the Ed-dimensional im- 
purity and along the "perpendicular" directions. The 
case Ed = is associated with point-like defects, and ex- 
tended parallel linear (planar) defects are related to the 
cases Ed = 1(2). Generalizing Ed to non- negative real 
numbers it may be interpreted as an effective fractal di- 
mension of a complex random defect system 6 ". The criti- 
cal behavior of such a model was examined by means of 
the renormalization group (RG) method 3 * 4 *^. A double 
expansion in both s = 4 — d, Ed was suggested and RG 
functions were calculated 3 - to order e, Ed ; qualitatively, 
the crossover to a new universality class in the presence 
of extended defects was found. These calculations were 
extended to the second order in Ref. 4 . Here, it was ar- 
gued, that the Harris criterion is modified in the presence 
of extended impurities: the randomness is relevant, if 



£ d > d . 



(2) 



For point-like disorder Ed = the Harris criterion (JIJ is 
restored from J2J. In particular, J2J) defines for each value 
of m a lower marginal defect dimension £™ ars , above 
which the critical exponents are influenced by disorder. 
Note that for a negative value of £™ arg any amount of 
quenched impurities will induce new critical behavior as 
far as Ed is always positive. 

Taking the best known estimates from a 6-loop RG ex- 
pansion for the exponent v v of the different m-component 
systems^ 5 ", one finds that the disorder with extended de- 
fects is relevant for d = 3 over a wider range of m than the 
point defect disorder. We show the lower marginal value 
£ marg |. nege systems in Fig. ^ Its asymptotic value for 



2 



large m is e™ ars (m = oo) = 1. This estimate is easily ob- 
tained using the value of the correlation length critical ex- 
ponent for d — 3 spherical modei-§ v (^l = 3 5 m — oo) = 1 
in ©. 
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FIG. 1: The marginal value of £™ arg for 3-dimensional m- 
component magnetic systems as function of m. Disorder is 
relevant in the region of Ed, m plane above the £™ ars curve. 
The solid line is obtained by substituting into the 6-loop 
results 15 for the critical exponents of pure w- vector magnet. 
The dashed line shows our present results, obtained in 2-loop 
approximation. 

Although the RG functions of Dorogovtsev's modefi 
were obtained in Rcf. 4 within the two loop accuracy, a nu- 
merical analysis was provided only to order e, Ed- For the 
class of extended-defect systems with cubic anisotropy, 
the second order E,£d RG analysis has been performed 7 . 

A related model with long-range-correlated quenched 
disorder, has been proposed in the work of Weinrib and 
Halperin^ and was recently studied in Refs>i2ii^. This 
model is characterized by a correlation function that has 
a power law decay g(r) ~ r~ b with distance r. This type 
of disorder has a direct interpretation for integer values 
of b. Namely, the case b = d corresponds to point-like 
defects, while b = d—l(b= d — 2) describes straight lines 
(planes) of impurities of random orientation, in contrast 
to the parallel extended impurities, considered above. 
The results of Weinrib and Halperin were confirmed by a 
Monte Carlo study 9 of the 3-dimensional random Ising 
model with linear defects of random orientation. An 
extension of this model describes d-dimensional systems 
with quenched disorder that is completely correlated in 
£d dimensions and long-range correlated in the remaining 
d — £d dimensions'^. This model applies when parallel 
dislocations in a crystal induce an elastic stress field in 
the direction perpendicular to the dislocation axis. The 
stress field may couple to the local critical temperature 
and play the role of a random potential. Real linear dis- 
locations are described by Ed = 1. The case Ed = 
reproduces the model of Weinrib and Halperin. The 
renormalization group treatment has been performed in 



a one-loop expansion within triple (e, Ed, a = d — b)H and 
(e, Ed,S = 4 — expansions. 

Although systems with extended quenched defects are 
subject to a number of studies, a complete classification 
of the critical behavior of magnets of this type has yet to 
be established. The results that have been obtained so far 
using second order double expansions are of qualitative 
character only. The series of the RG functions are known 
to be asymptotic at best, and to obtain reliable quantita- 
tive results, appropriate resummation techniques should 
be applied. Such have proved to be fruitful for the anal- 
ysis of point-like quenched disorder^S^. In our previ- 
ous paper— we have already applied this analysis to the 
3-dimensional magnetic systems with extended defects. 
Previously, resummation techniques were not exploited 
for these systems, apart from the special case d = 3 and 
6d = 1, corresponding to linear defects^. There, critical 
exponents were calculated to second order in an expan- 
sion in s, Ed and Pade-like approximants were used to 
provide numerical estimations. In our present paper we 
explore a wider region of the phase diagram shown in 
Fig. n an d apply resummation schemes to estimate nu- 
merically the critical exponents of extended-defect sys- 
tems both in the case of space dimension d — 3 and in 
general dimensions, using the two loop RG functions of 
RefsM. 

Our paper is organized as follows. In the following 
section II we present the model. The field-theoretical 
analysis of critical behavior is reviewed in section III, and 
expressions are given for the RG functions. Resummation 
is applied to the RG functions both at d = 3 and in the 
case of general dimensions and numerical estimates for 
the critical exponents are presented in section IV followed 
by conclusions in section V. 

II. THE MODEL 

We consider the model^^ of an m-vector magnet 
with e^-dimensional impurities, each extending through- 
out the system along the coordinate directions symbol- 
ized as afii, whereas in the remaining d — Ed dimensions 
they are randomly distributed. The effective Hamilto- 
nian of the model reads^: 

H = J d d x[^ 2 + V(x))$ 2 (x) + (Vj(x)) 2 

+ ao (V||^)) 2 + ^(^(x)) 2 ]. (3) 

Here, <f> is an m-component vector field: <f> = {0 1 • ■ • </>™ 1 }, 
/xo and Uq are the bare mass and the coupling of the 
magnetic model, a is the bare anisotropy constant, V|| 
and stand for differentiation in the coordinates x\\ 
(along Ed directions) and x± (along (d — Ed) directions) 
respectively, V(x) represents the impurity potential. The 
probability distribution for the impurities is defined to 
yield: 

((V(x))) = 0, 



3 



((V(x)V(y))) 



(4) 



Here, ((...)) stands for the average over the potential dis- 
tribution, (-i>o ) is a positive constant proportional to both 
the concentration of impurities and the strength of their 
potential. The extended defects make the space coordi- 
nate anisotropic; the anisotropy is parameterized by the 
constant clq. It can initially be assigned a value of unity, 
but will be renormalized away from this value in the scal- 
ing limil£. 

As far as the disorder is considered to be quenched, 
the free energy of the system is obtained as follows^: 



T=-kT({hxZ c )), 



(5) 



where Z c is the partition function for a given (quenched) 
configuration of impurities, and the average is performed 
over the impurity probability distribution Q). To avoid 
averaging the logarithm of the partition function in JSJ), 
the standard procedure is to use the replica trick^i . This 
amounts to formally implementing the identity In Z = 
lim n _ >o{Z n — l)/n by analytic continuation in n. After 
performing the average the n-replicated effective Hamil- 
tonian then reads& 



n . 

= E / dM^ote + (Vx<Mz)) 2 
a =\ J z 

+ a (V||to) 2 ] + ^(to) 2 } 



(6) 



vq 
2 



]T Jd d xJ d d y6 d -^(x x -y x )^ a (x)^(y). 

OL.J3— 1 

Here, Greek indices denote replicas. The last term intro- 
duces an additional coupling vq < 0. It is present only 
for non-zero dilution and is directly responsible for the 
effective interaction between replicas due to the presence 
of impurities. The replica limit n — > is implied for all 
quantities calculated in the following. 

The model described by the effective Hamiltonian © 
has a rich scaling behavior. An introduction to its criti- 
cal exponents and a derivation of their scaling relations 
may be found in Refs^^. Here, we summarize the main 
quantities that describe the critical behavior of the model 
© . Due to the spatial anisotropy two correlation lengths 
exist, one perpendicular and one parallel to the extended 
impurities direction: and £| | ) . As the critical tem- 
perature T c is approached, their divergences are charac- 
terized by corresponding critical exponents v±_ , Vu: 



1*1" 



1*1" 



(7) 



where t is the reduced distance to the critical temperature 
t = (T — T c )/T c . The correlation of the order parameter 
fluctuations in two different points depends on the ori- 
entation of their distance vectoA Thus, the critical ex- 
ponents rj± and n | , that characterize the behavior of the 
correlation function in the directions, perpendicular and 
parallel to the extended defects, must be distinguished. 



On the other hand, as far as the interaction of all order 
parameter components with defects is the same, the sys- 
tem susceptibility is isotropic and can be expressed by 
the pair correlation function^: 



x(ki,*n,t) = ir 7 ff(- 



fc ll 1*1 1 

i K n 



(8) 



In (jSJ), fen, k± are the components of the momenta along 
Ed and d — Ed directions, respectively, 7 is the magnetic 
susceptibility critical exponent and g is the scaling func- 
tion. For the above introduced critical exponents the 
following scaling relations hold^*&: 



7 = (2- r)±)i>± = (2 - r?||)^||. 



(9) 



The critical exponent a of the specific heat is related to 
z/j_j v\\ by another scaling relation that differs from the 
ordinary one^: 



a = 2 - (d - e d )vx - E d v\\. 



(10) 



All the other scaling relations are of the standard form 22 . 
This implies that one should calculate at least three inde- 
pendent exponents (e.g., M 1 , fj_ f 7) instead of two, as in 
the standard case, to find the others by scaling relations. 



III. RENORMALIZATION GROUP FUNCTIONS 

To describe the long-distance properties of the model 
© near the second order phase transition we use the 
field-theoretical renormalization group (RG) method 23 . 
In this approach one observes the behavior of the system 
under a scaling transformation that rescales the length 
e.g. by a factor I. At the critical point the system as- 
sumes scale invariance. Away from the critical point the 
change of couplings uq, vq — > u, v and the anisotropy con- 
stant clq — > a under rescaling defines a flow in parametric 
space, expressed by RG functions: 



A, 



du 1 



a - dv 1 



C^l.. (H) 



Here, £ is the rescaling factor and the notation |o indicates 
differentiation at fixed bare parameters. The bare field 
<f> and the bare mass /iq are related to the renormalized 
field if and mass fj, by: 



2 ry 2 



(12) 



The fixed points (FPs) u*,v* of the RG transfor- 
mation are the solutions of the system of FP equa- 
tions: /3 u (u*,v*) = 0, p v (u*,v*) = 0. A FP is sta- 
ble if the eigenvalues of the stability matrix Bij = 
dp Ui (u* ,v*)/df3 Uj (u*,v*) have positive real parts. The 



4 



stable FP that can be reached starting from the initial 
values of the coupling constants (in our case, u > 0, 
v < 0), corresponds to the critical point of the system. 
At this point, the magnetic susceptibility and correlation 
length critical exponents are given by by the relations: 



7 



1 = 1- 



70 2 



2-70 
'I 1 = 2-7 2 - 70 , 

v\\ = (l - yW> 



(13) 



while the expressions for the RG functions 7^,7^2 read: 



70 



dial 



14> 2 



dlnZ-2 1 
dint 1 



■70- 



(14) 



The other critical exponents can be obtained from the 
scaling relations ©, ijTU|) . 

In order to derive the quantitative characteristics of 
the critical behavior of magnetic systems with extended 
impurities, we analyze the 2-loop RG functions, obtained 
in* in the minimal subtraction scheme. Unfortunately, 
that paper contains some numerical errors, as pointed 
out in Ref£. The appropriately corrected functions read: 



Pu/U 



—e - 
1 

12* 



(m + 8) n (3m + 14) 2 
-u -2v- - -u 



6 



12 



|(llm + 58) + (m-4) ^ 
3 3(e + £ d ) 



144 



328 + 32- 



0v/v = 

+ 

70 = 
70 2 



-e - e d - -v 
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£d 
to + 2 



18 



3 



-u - -v 
6 



m + 2 



5 m + 2 o 

-« 2 ; (15) 



12 3 



36 36 
to + 2 1 



72 

to + 2 



6 

to + 2 

\ 

24 



v 

3 



u 2 - 24v 2 



Ed 



£ + £d 



(16) 
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D 1; H (2 + m)vu. 
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In the case of Ed = the dimension enters the minimally 
subtracted RG /3-functions l|15|) only in the form of the 
engineering dimension e. For fixed dimension e = 1 the 
flow Ijllfl with the /3-functions l|15|l can then be directly 
evaluated 2 ^ 2 ^. Also for non-zero £d no singularities are 
introduced by the dimension dependence. We thus pro- 
pose to extend the approach of direct evaluation to the 
RG functions of the present model, i.e. to treat them 
directly at d — 3 (e = 1) for different fixed values of the 
(fractal) defect dimensionality Ed- 

The expansions of the RG functions in powers of the 
coupling are known for a number of models with high 



accuracy. However, the very idea of perturbation the- 
ory is that the result may successfully be approximated 
by accounting for higher order contributions. In the 
field-theoretical RG approach the series appear to be 
divergent; moreover, they are characterized by a facto- 
rial growth of the coefficients implying a zero radius of 
convergence^. A simple method to manage this problem 
is to truncate the series optimally, i.e., to account only for 
those first several terms that do not show the divergence. 
To take into account the higher order contributions, the 
application of special tools of resummation is required 26 . 

There remains the principal question about the Borel 
summability of the perturbation series for a given model; 
at present such a proof has been given only for the 4> 4 
theory with one coupling 2 ^. The field-theoretical RG se- 
ries for models with several couplings were analyzed as 
if they are asymptotically divergent without a proof of 
this property. Moreover, there exists strong evidence of 
possible Borel non-summability of the series obtained for 
disordered models 2 ^. 

Recently, it was demonstrated analytically, that the 
RG functions of the d = random Ising model are 
Borel summable, if a special resummation technique is 
exploited 2 ^, the application of this technique to the d = 3 
RG series of the random Ising model allows to restore the 
convergence of the results obtained in the massive RG 
scheme^S and in the minimal subtraction scheme^!. 

In our case, the question about the summability of the 
series I|15f) - (ll6t remains open. Nevertheless, we apply 
various kinds of resummation techniques to obtain re- 
liable quantitative results. In the following section, we 
describe the procedure of the Chisholm-Borel resumma- 
tion technique, which proves to be most effective in our 
case, and we present numerical results for the critical 
behavior of systems with extended defects both in the 3- 
dimensional case and in the case of general (non-integer) 
dimension. 



IV. THE RESULTS 
A. Resummation procedure 

The two-variable Chisholm-Borel resummation 
techniqueiii^ consists of several steps. To explain, how 
it is applied to the RG functions let us denote any of the 
expressions <|15H <|16[) by / = f(u, v). First, we construct 
the Borel image of the initial RG function /: 



/=£< 



Ea 
1 



r(*+i + i) ' 



where T(x) is Euler's gamma function. Then, the Borel 
image is extrapolated by a rational Chisholm 32 approxi- 
mant [K/L](u,v). One constructs this ratio of two poly- 
nomials of order K and L such that its truncated Taylor 
expansion is equal to that of the Borel image of the func- 
tion /. The resummed function is then calculated by an 
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inverse Borel transform of this approximant: 



r 



dtexp(-t)[K/L](ut,vt). 



(17) 



There is a lot of possibilities to choose a Chisholm ap- 
proximant in two variables. The most natural way is to 
construct it such that, if any of u or v is equal to zero, the 
familiar results are obtained for reduced model. Here, for 
the Borel images of the /3-functions we have chosen the 
following approximant with a linear denominator: 







chis 



bo.o + bi.out + bo.ivt + b lt iuvt 2 

1 + Ci nUt + CQiVt 



(18) 



Note, that the polynomial in the numerator is chosen to 
be symmetric in the variables u and v. As far as the 
Chisholm approximants enter the integrals (|17fl . the ex- 
pression under integration may contain poles; in this case 
an analytic continuation is needed, taking the principal 
values of the integral. Below we display results only for 
those situations that do not require such an analytic con- 
tinuation. 

To calculate the values of the critical exponents v±,U\\, 
7 we substitute ((Ttjj) into |TSJ), apply the resummation 
procedure to the resulting series, and, finally, substitute 
the obtained values of stable fixed point coordinates in 
the above expressions for the critical exponents. Note, 
that while we use the Chisholm approximant in the form 
(|18fl to analyze the series for u±,'y, the symmetry of the 
series for allows us to use the Chisholm approximant 
in the form: 



1 + C\fiUt + CQivt 



(19) 



In the following we apply the above resummation tech- 
nique to analyze the critical behavior of the model © 
both at space dimensionality d = 3 and in general di- 
mensions. 



B. 3-dimensional systems with extended impurities 

Let us first consider the model Eq. JJjJ for the three- 
dimensional case. Fixing the value e = 1 (i.e. d = 3) 
and treating Sd as a variable parameter, we look for the 
common zeros of the resummed functions f3 u /u and f3 v /v. 
Let us note that there are two possible scenarios for the 
change of FP stability under disorder introduced to the 
system. This is due to the fact that the universality class 
of the weakly diluted m-vector model with point-like de- 
fects {ed = 0) depends in a crucial way on the order 
parameter dimensionality m. There exists a marginal 
value m c such that for m > m c the critical exponents 
of the m-vector model remain unchanged by point-like 
defects whereas for m < m c they split to new values. 
Discussing the Harris criterion we mentioned, that for 
integer m point-like defects lead to a new universality 



class only for the Ising model (m = 1). The present esti- 
mates for m c based on the resummation of six-loop RG 
functions definitely imply m c < 2: m c = 1.942 ± 0.026 33 
and m c = 1.912 ± 0.00421 

For m < m c , at Ed — the random FP (u ^ 0, v ^ 0) 
is stable and physically accessible. Any small increase of 
the parameter £d leads to a shift of the stable FP value 
and therefore to new critical exponents, compared to the 
pure ones (see Fig. [21 a)- For m > m c the situation is 
different. Here, the point-like (Ed — 0) disorder is not rel- 
evant, and the pure fixed point (u* ^ 0, v * — 0) is stable, 
whereas the random fixed point lies in the unphysical re- 
gion and is unstable. Switching on the parameter Ed leads 
to a shift of the random FP, and at some value of this 
parameter the random FP passes to the physically ac- 
cessible region, at the same time it interchanges stability 
with the pure FP, i.e. it becomes stable (see Fig. |3 b). 
This corresponds to a crossover to another universality 
class with different critical exponents. 



, s d = 




a) m < m c 



b) m > m c 



FIG. 2: The two qualitative scenarios for the change of fixed 
point stability. Circles mark unstable fixed points, squares are 
the stable ones. The polymer FP (u* — 0, v* ^ 0) is stable 
but unphysical. a) The spin dimensionality of the magnetic 
system is below the marginal one. The random FP (u* ^ 
0, v* ^ 0) is stable at Ed = 0. As the parameter Ed increases, 
the mixed fixed point remains stable and moves within the 
physically accessible region, b) The spin dimensionality of 
the magnetic system is above the marginal one. The pure FP 
(u* 7^ 0, v* — 0) is stable at Ed = 0, whereas the random FP 
is unstable and lies in the unphysical region. As Ed increases, 
the mixed fixed point moves towards the physical region and 



at £ d 



becomes stable. 



The the stable FP coordinates obtained for the 3- 
dimensional m-component magnetic systems with m 
l.2..'i. -". are presented in Table [H One sees that for m = 
2,3,4 (i.e. for m > m c ) the pure FP [u* ^ 0,v* = 0) 
looses its stability for some marginal value £™ arg . As 
mentioned in the introduction, £™ alg can be calculated 
from the correlation length critical exponent of the pure 
magnet with help of the generalized Harris criterion J5J. 
We will return to these estimates below. 

The results obtained for the resummed values of the 
critical exponents are presented in Table [n] We may 
draw some conclusions from these results. As said above, 
the case Ed = describes point-like quenched disorder, so 
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it reproduces the well known results^; in this case there 
are no longitudinal components of the critical exponents. 
For Ed 0, the relation 1/11 > v± holds for every Sd and to. 
This may be explained in the following way: the extended 
defects cut interacting paths of spins perpendicular to the 
extended-defect direction, so in the parallel direction the 
fluctuations are stronger and the correlation length more 
sharply diverges. 
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FIG. 3: The critical exponent of the correlation length per- 
pendicular to the ed-dimensional "lines" of extended defects 
as function of the parameter Ed for three-dimensional m- 
component spin systems. 
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FIG. 4: The critical exponent of the magnetic susceptibility 
of three-dimensional m-component spin systems as function 
of the parameter Ed- 

In Figs. |3lE]we plot the values of the critical exponents 
v±_ and 7 for TO-component systems with m = 1,2,3,4 
as functions of the parameter Ed- They clearly represent 
our discussion above: as we can see, in the case of the 
Ising model (to = 1), any small increase of Ed results in a 
change of the critical exponents, whereas for m = 2, 3, 4 
the critical exponents remain constant and equal to the 



corresponding exponents of pure model, until Ed reaches 
its marginal value, and only for Ed > £™ arg the values 
of exponents increase. Note, that as the parameter Ed 
changes from to 1, the resulting change in the values of 
the critical exponents are sufficiently large to be experi- 
mentally observed (e.g., the value of v± for Ising model 
increases by 13 %). Thus, the influence of extended de- 
fects on the critical behavior of three-dimensional mag- 
netic systems can be observed experimentally. 

Let us return to the estimates of the marginal value 
£™ arg , at which the crossover to a new universality class 
occurs. In Fig. ^ we plot the estimates for £™ arg 
as obtained from six-loop results^ for the correlation 
length critical exponent v p (m) of the pure to- vector mag- 
net: e^ s (m = 1) = -0.173; e™ IS (m = 2) = 0.016; 
e™ r S(m = 3) = 0.172; £™ arg (m = 4) = 0.300. The RG 
functions of the m-vector magnet with extended impu- 
rities exploited in this paper (expressions l|15 |) -l|16 p ) are 
currently known with two-loop accuracy. In this approx- 
imation we get for the marginal value e™ arg (TO = 1) = 
-0.105; e^ arg (TO = 2) = 0.077; e™ arg (m = 3) = 0.222; 
e? arg (m = 4) = 0.337. The phase boundary of the crit- 
ical behavior of 3-dimcnsional m-component magnetic 
systems in the m, e^-plane is shown in Fig. ^ by the 
dashed curve for the two-loop approximation. It is re- 
markable that the resummed two-loop and six-loop re- 
sults for £™ arg lead to a very similar phase diagram. In 
contrast to this the first order e, e^-expansion predicts 
a qualitatively different behavior—: the random FP is in 
the physical region for e > Ed, and randomness is rele- 
vant for to < to c = 4(e + 2£d) /{s — £d), while for to > m c 
the random FP lies in the unphysical region and the pure 
FP is stable, so for e > Ed the stability changes at m c . 
For e < Ed the only physical fixed point is random and is 
always stable. From these first order results one would 
conclude, that the disorder is relevant for every to < 4 
and positive Ed- 

Another interesting question concerns the existence of 
an upper marginal value for the defect dimensionality 
Ed- This question has not been raised in previous works, 
where the double e, e^-expansion has been exploited. In 
our analysis, we observe the disappearance of a stable 
reachable FP for Ed slightly above 1. This, in principle, 
may serve as evidence of an upper boundary for Ed above 
which the second order phase transition ceases to exist. 
However, this behavior was observed in the region of cou- 
plings u, v where the Chisholm approximants (|18fl that 
enter the expressions for the resummed /3-functions l|17|) 
have poles on the positive real axis. This excludes a def- 
inite answer about the presence and stability of the fixed 
points for high Ed- On the other hand, simple reasons 
support the idea of an upper boundary for Ed- Inter- 
preting Ed as the fractal dimensionality of the defects, 
it is clear, that it cannot exceed the dimensionality of 
the embedding space, d = 3. Furthermore, one expects 
physically, that extended defects of large dimension (e.g. 
parallel planar defects with Sd = 2), that extend through- 
out the system, divide the system into non-interacting 
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regions and thus inhibit ferromagnetic order. 



C. Systems with extended defects in general 
dimensions 

Let us now consider the case of non-integer space di- 
mension d. The concept of non-integer space dimension 
is common in the theory of critical phenomena. Consid- 
ering d as a continuous variable is not purely formal; for 
the critical behavior of spin systems placed on the sites 
of self-similar fractal lattices^ it has a direct geometrical 
interpretation. 

Considerable effort has been made to calculate the crit- 
ical exponents for the second order phase transition for 
general (non-integer) values of d (see Refsi^SL^ for re- 
view). Most papers study the Ising model. The criti- 
cal exponents of the pure m-vector model in general di- 
mensions have been calculated in the field-theoretical RG 
scheme for the case m = 1 (Ising-like systems) 3 ^! and 
for m = 2, 3, 42£. Numerical results for the weakly di- 
luted Ising model are presented in Refi^. 

To analyze the critical properties of the model © for 
general dimensions, the 3d approach, used above for the 
RG functions of the present model is generalized. We 
treat the RG functions directly at different fixed values 
of (non-integer) space dimension d and defect dimension- 
ality £ d . 

Again we are looking for the marginal value of the pa- 
rameter Ed for models at different non-integer space di- 
mensions. To make use of the generalized Harris cri- 
terion in this case, we take the known estimates 
from a 3-loop RG analysis of the critical exponent 
v p {d) for the pure m-vector model at m = 1— and 
rn = 2,3,4^. In accordance with the Mermin-Wagner- 
Hohenberg theorem^, the lower critical dimension of the 
m-vector model is d^ = 2 for m > 2, while in the case of 
the Ising model (m = 1) one has d^ = l 22 (for d < d^, 
there is no spontaneous magnetic ordering at T > 0). 
In order to take into account these conditions, and also 
to reproduce the exact Onsager valued for the 2d Ising 
model (y = 1), a more sophisticated resummation pro- 
cedure was exploited in Refsi^LiS, and we will use these 
results. 

Substituting the known data™ i2£ for v(d) into the gen- 
eralized Harris criterion (J2J we find the marginal value 
of the parameter Ed for different m-component systems 
as function of space dimension, see Fig. \5\ This figure 
serves as a phase diagram in parameter space id, Ed)- It 
should be noted, that, as follows directly from the Harris 
criterion (JTJ, one has £™ alg = at d = 4 (the upper crit- 
ical dimension, where the Gaussian FP (u* = 0, v* = 0) 
is stable and v = 1/2), whereas at d = dz,, one has 
£™ arg = di,, independent of m. Note, that for m = 1, 
2 there exists a range of space dimensions d, at which 
e™ arg < 0, and, thus, the presence of any extended de- 
fects causes changes to the critical behavior. In the region 
marked as I in Fig. [S] disorder with extended defects is 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 



FIG. 5: The phase diagram of m-component magnetic sys- 
tems with extended defects in the d, £,j-plane. The lower 
marginal value of e d for m = 1 (squares), 2 (circles), 3 (tri- 
angles), 4 (diamonds) is shown. In region I disorder with 
extended defects is irrelevant for all m > 1; region II: disor- 
der is relevant only for the Ising model (m = 1); region III: 
disorder with extended defects is relevant for systems with 
m = 2,3, 4. The solid diagonal line corresponds to Ed = d: 
the region IV corresponds to non-physical situations. 



irrelevant for all the systems under consideration; in re- 
gion II it becomes relevant for the Ising model, passing 
from the region II to III it gradually becomes relevant 
for the models with m = 2, 3, 4. As discussed above, the 
concept of an upper marginal value for the parameter Ed 
naturally arises. The defect dimensionality Ed cannot ex- 
ceed the dimension of the embedding space d, and thus 
cannot take values from the region marked as IV in Fig. 
□ 

To get the values of the critical exponents of the set of 
systems with extended defects at general dimensions d, 
we exploit the same 2-variable Chisholm-Borel resumma- 
tion technique as for the 3-dimensional case above. But 
now we calculate the RG functions l|15|) - (|16|) for fixed 
non-integer e = 4 — d. We present here only the results 
obtained for the extended-defect Ising model (m = 1), as 
far as it is most sensitive to introducing defects in a wide 
region of space dimensionality d, as can be observed from 
Fig. E 

Before discussing the results obtained for general d let 
us note, that their accuracy depends on the distance e = 
4 — d from the upper critical dimension d = 4^^^^- 
for larger expansion parameters e also the fixed point 
values of the coupling constants grow leading to expan- 
sions which are difficult to estimate, even when resumma- 
tion procedures are applied. We demonstrate this for the 
d — 2 and d — 3 Ising model with point- like quenched dis- 
order {Ed = 0). The Onsager solution of the d = 2 pure 
Ising model 41 brings about a logarithmic divergence of 
the specific heat, i.e. an exponent a p = 0. With the 
Harris criterion Q the well known irrelevance of this 
type of disorder for the d = 2 Ising model follows. The 
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FIG. 6: The critical exponent of the correlation length per- 
pendicular to the Ed-dimensional "lines" of extended defects, 
as function of dimension d for Ising systems at values of the 
parameter Ed = 0, 0.2, 0.4, 0.6, 0.8, 1 from below. 



y 1.6 - 




2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 

of 



FIG. 7: The critical exponent of the magnetic susceptibility 
as a function of the dimension d for Ising systems at values 
of the parameter Ed = 0, 0.2, 0.4, 0.6, 0.8, 1 from below. 



correlation length critical exponent in this case equalsai 
v p = 1, leading to £™ arg = by applying the generalized 
Harris criterion (J2J. However, in our 2-loop calculations 
refined by the Chisholm-Borcl resummation technique we 
find for the d = 2 Ising model correlation length and 
magnetic susceptibility critical exponents v p = 1.162 and 
7 P = 2.042 correspondingly. These values differ by 16 % 
from their exact values 41 ^ v = 1, j p = 7/4. This implies 
a marginal value of the extended defects dimensionality 
of £™ arg = 0.278. For d — 3, we reproduce the two- 
loop results of the random Ising model critical exponents, 
obtained in the minimal subtraction scheme in Rcf. 25 
Ux = 0.665, 7 = 1.308. They differ from the most accu- 
rate six-loop estimates^ v± = 0.678(10), 7 = 1.330(17) 
by only of 2 %. Again, for the marginal value of the de- 
fect dimensionality we find £™ alg = —0.105, the six-loop 



value being £™ alg = —0.173 (see Fig. (1)). To conclude, 
we expect the accuracy of our results for the critical expo- 
nents to be of the order of several percents for d = 3 and 
to increase to the order of ~ 10% — 20% for d = 2, (recall 
that the approach gives the exact data for d = 4). For 
this reason we do not show any results for d < 2, taking 
that the confidence interval for them is even larger. On 
the other hand, the two-loop approximation restricts our 
investigation to dimensions d smaller than 3.4; namely, 
for d > 3.4 the rational approximants that enter the inte- 
grals for the resummed functions have poles that exclude 
reliable calculations in this case. 

The values of the critical exponents v±_ and 7 are plot- 
ted as functions of d parameterized by Ed in Figs. El 
|7| . The parameter Ed of the different curves changes 
from to 1 by steps of 0.2. For values of d below 2.1 
the critical exponents remain equal to those of the pure 
model for small Ed, until Ed reaches its marginal value (of 
about Ed — 0.27 in the present two-loop calculations); 
then, a crossover to the diluted regime takes place, and 
the exponents increase continuously with increasing the 
parameter Ed at fixed space dimension d. We interpret 
this by noting that an increase of Ed leads to an effective 
decrease of the space dimension implying a sharper diver- 
gency of the correlation length (characterized by a larger 
value of v). The results for z/11, not shown here, display 
the same tendencies under changes of the parameters d, 
Ed- Let us note however, that the relation vu > v± holds 
for all d, Ed, as observed for d = 3 and explained in the 
previous subsection. 



V. CONCLUSIONS 

In our study we provide numerical estimates for the 
critical exponents of m-component magnetic systems in 
the presence of extended macroscopic defects. The impu- 
rities are envisaged as Erf-dimensional objects, extending 
throughout the system, whereas in the remaining d — Ed 
dimensions they are randomly distributed. The presence 
of these impurity "lines" introduces anisotropy into the 
systems; as a result, the parallel and transverse corre- 
lation lengths naturally arise and are described in the 
vicinity of the critical point by the respective correlation 
length critical exponents^. Such systems are not covered 
by the original Harris criterion 4 ; disorder with extended 
impurities is relevant over a wider range of m and d than 
the point defect disorder (see Figs. Q][5J|. 

Although the model of m-component magnetic system 
with extended quenched defects has attracted much at- 
tention and serves as a subject of study in a number 
of works, most results are based on the double expan- 
sion in parameters s,£d and are rather of a qualitative 
character. We analyze the RG functions of the model, 
obtained in the minimal subtraction scheme^, treating 
them directly for fixed d and fixed parameter Ed- To 
evaluate the RG functions we apply appropriate resum- 
mation techniques, that have proven to be fruitful in the 
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analysis of point-like quenched disorder—. In our case, 
the question about the summability of the corresponding 
RG functions is open; nevertheless, in order to obtain reli- 
able quantitative results, we apply a simple two-variable 
Chisholm-Borel resummation technique. The values of 
the stable fixed point coordinates and critical exponents 
for 3-dimensional systems are extracted, as well as the 
estimate for the lower marginal value for parameter Ed is 
obtained. The case e c i = describes point-like quenched 
disorder and reproduces well-known results^. For > 
it was found, that the relation v\\ > u± holds for every 
£d and to, which describes the following physical situa- 
tion: the extended defects cut interacting paths of spins 
perpendicular to the extended-defect direction, so in the 
parallel direction the fluctuations are stronger and the 
correlation length more sharply diverges. The only nu- 
merical results for the exponents of the model © were 
obtained so far for Ed — 1 using the Pade-analysis^, they 
are shown in the last line of Table ITT1 We note however, 
that our present resummation technique has shown its ef- 
ficiency and accuracy in studies of models with point-like 
structural disorder— and is generally known to provide 
more reliable estimates as compared to the simple Pade- 
analysis. 

The critical behavior of the 3-dimensional Ising model 
(m = 1) is changed by an extended impurity character- 
ized by any Ed ^ 0, whereas for to = 2, 3, 4 the extended 
defect disorder is relevant only for Ed > £™ arg - We es- 
timate the marginal value of the parameter Ed also for 
magnetic systems in general (non-integer) space dimen- 
sions. In Figs. EJ{7|we plot the results for the Ising model 



in general dimensions, which appears to be the most "sen- 
sitive" to introducing extended defects in a rather wide 
range of space dimensionality d. 

The concept of an upper marginal value of Ed natu- 
rally arises in studies of magnetic systems with extended 
impurities. Interpreting Ed as a fractal dimensionality 
of defects, it is clear, that it cannot exceed the dimen- 
sionality of embedding space, d. Moreover, physically 
one may expect, that extended defects of large dimension 
(e.g. parallel planar defects with Ed = 2), each extend- 
ing throughout the system, will divide the system into 
non-interacting regions and thus prevent it from ferro- 
magnetic ordering. 

Unfortunately, there is no simulational investigation 
of systems with parallel extended defects that we know 
of 42 . However, our calculations bring about an essential 
change in critical exponents due to presence of extended 
impurities (e.g. vy_ for the Ising model with extended 
impurities increases by 13 % as the parameter Ed changes 
from to 1). We hope that the influence of this kind 
of disorder on the critical behavior of three-dimensional 
magnetic systems may create some interest for numerical 
simulations or experimental measurements. 
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TABLE I: The coordinates of the stable fixed points of the 
RG equations for 3-dimensional magnetic systems with ex- 
tended Ed-dimensional impurities for different m. Note, that 
the pure fixed point (u* 7^ 0, v* — 0) looses its stability for 
some marginal value e™ arg (m). 
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TABLE II: The critical exponents for 3-dimensional magnetic 
systems with extended e-dimensional impurities for different 
m. For comparison the last line shows the results of Refi—. 



